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ABSTRACT 

We discuss the formation of pulsars with massive companions in eccentric orbits. We 
demonstrate that the probabihty for a non-recycled radio pulsar to have a white dwarf 
as a companion is comparable to that of having an old neutron star as a companion. 
Special emphasis is given to PSR B1820-11 and PSR B2303+46. Based on population 
synthesis calculations we argue that PSR B1820-11 and PSR B2303+46 could very 
well be accompanied by white dwarfs with mass ^ 1.1 M0. For PSR B1820-11, how- 
ever, we can not exclude the possibility that its companion is a main-sequence star 
with a mass between ^ 0.7 M© and ^ 5Mq. 
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1 INTRODUCTION 

High mass binary pulsars are binaries in which a radio pul- 
sar is accompanied by an unseen companion that is approxi- 
mately as massive as the pulsar. Although the companion is 
generally believed to be another neutron star, it appears to 
be likely that some of the high-mass binary pulsars are not 
accompanied by neutron stars but are instead accompanied 
by white dwarfs or low-mass main-sequence stars. 

PSR B1820-11 (J1823-1115), for example, is a radio pul- 
sar in a binary with an orbital period Porb ~ 357.8 days and 
an eccentricity e « 0.795 (Lyne & McKenna 1989). The de- 
rived mass function / — 0.068 Mq indicates that the mass 
of the companion exceeds ~ 0.7 M©. The small spin-down 
age of the pulsar (r = 6.5 Myr) and its strong magnetic field 
(B ~ 6 X 10^^ G) suggest that it is a young - non-recycled - 
pulsar (Phinney & Verbunt 1991). 

Based on the mass function and the orbital eccentricity, 
Lyne & McKenna (1989) suggest that the most likely com- 
panion to PSR B1820-11 is a neutron star. The companion 
is unlikely to be a massive hydrogen star or a black hole 
due to the low mass function. Lyne & McKenna argue that 
the observed orbital eccentricity excludes the companion as 
a white dwarf or a low-mass helium star. Phinney & Ver- 
bunt (1991) suggest that the companion of PSR B1820-11 
may be a ~ 1 Mq main-sequence star, which will make the 
system a precursor of low-mass X-ray binary. 



Also PSR B2303-I-46 may contain a non-recycled pul- 
sar (P = 1.06 s, r = 30 Myr, B = 8 x 10^^ G) in a binary 
with Porb = 12.34 day and e fa 0.66 (Stokes et al. 1985). 
The total mass of the binary is 2.64 ± 0.05 M© (Thorsett & 
Chakrabarty 1999) and no optical counterpart was detected 
down to P = 26™ (Kulkarni 1988). Recently van Kerkwijk 
& Kulkarni (1999), however, discovered an object which co- 
incides with the timing position of PSR B2303-I-46 and the 
properties of this object are consistent with those of a mas- 
sive white dwarf with a cooling age close to the age of pulsar. 

With population synthesis calculations we demonstrate 
that the birth rate of binaries in which the neutron star is 
born after the white dwarf is comparable to the birth rate 
of binaries with two neutron stars. In the former case the 
binary will have a high eccentricity eccentric and the radio 
pulsar will be accompanied by a white dwarf with a mass 
> 1.1 Mq. In this paper we argue that both PSR B1820-11 
and PSR B2303-I-46 systems may be formed via such a sce- 
nario (see §2.2). The calculations, however, do not rule out 
that the companion of PSR B 1820-11 is either another neu- 
tron star or a main-sequence star. In the latter case, the 
mass of the main-sequence star would most likely be be- 
tween ~ 0.7 M© and ~ 5Mo. 
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2 THE BINARY COMPANIONS OF 
PSR B1820-11 AND PSR B2303+46 
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2.1 Hiding a main-sequence star 

At a distance of 6.3 kpc and a height above the Galac- 
tic plane of 100 pc (Taylor, Manchester, & Lyne 1993), 
PSR B1820-11 is heavily obscured by the interstellar 
medium. A star with an absolute magnitude of My ~ 
— 4!"0 could easily be missed in a survey with a limiting 
magnitude of 20!"0 (assuming an interstellar absorption of 
Av = 1!"6 kpc"^). In the USNO-A vl.O catalog (Monet et al. 
1997), based on the Palomar Sky Survey, the object closest 
to PSR B1820-11 is located at a distance of 5.9 seconds of 
arc and therefore cannot be associated with PSR B1820-11. 
There is no indication for other forms of radiation (X-rays, 
infrared, etc.) from the direction of PSR B1820-11. These 
suggest that the companion can not be a very massive star 
or have a strong stellar wind. 

The presence of radio emission indicates that the rate of 
mass loss by the companion star is modest. An upper limit 
to M in the stellar wind of the companion can be estimated 
from the expression for the optical depth of the stellar wind 
for free- free absorption {e.g. , Eq. 16 in lUarionov & Sunyaev 
1975). By assuming that the wind of the accompanying star 
is transparent at A = 75 cm (and a temperature of the stellar 
plasma of 10* K) we obtain 



M< 4.810-^^^^(1 
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where Mtot is the mass of the binary. The estimate for M 
depends on the eccentricity; substitution of (1 -I- e) for (1 — 
e) in Eq. (^ gives the lower limit for the mass loss rate 
given that the radio pulsar is only visible at apocenter. We 
illustrate this limit in Fig. 1 for two typical combinations of 
Mtot and M. 

A binary with more or less similar orbital charac- 
teristics is PSR B1259-63 (Porb = 1237 ± 24 days, e = 
0.8699, Adopt = 10 ± 3Mq, Af = 5 x 10"* Mq yr^^ John- 
ston et al. 1992, 1994) which contains a radio pulsar and 
a Be-type star. The pulsar is visible at apastron but not 
at periastron. This is understood from Eq. (0) as an effect 
of shielding of the radio signal from the pulsar by the stel- 
lar wind. Another known radio pulsar in a long-period bi- 
nary is PSR J0045-7319 (Porb = 147.8 days, e = 0.808) 
with a Bl V (M = 8.8 ± 1.8 M©) companion with M < 
3.4x 10-"(uinf/wosc) Mq yr"^ (McConnell et al. 1991; Bell 
et al. 1995; Kaspi, Tauris & Manchester 1996). This pulsar 
is visible along its entire orbit due to a more tenuous wind of 
the companion star, which is typical for a lower metallicity 
star in the Large Magellanic Cloud. 

The absence of eclipses and of X-rays from accretion or 
shocks generated by the interaction between the pulsar wind 
with wind of the Be star (as are observed in PSR B1259-63, 
Grove et al. 1995) also indicates that the companion of 
PSR B1820-11 is not likely to be a massive star. An up- 
per limit to the mass-loss rate of the companion star in 
PSR B1820-11 may be set by measuring variations in the 
dispersion measure of the radio pulsar when it is close to 
periastron (see Melatos, Johnston, & Melrose 1995). 



2.2 Can the companion to a radio pulsar be a 
massive white dwarf? 

The eccentricity of a binary orbit can either be primordial or 
induced by the sudden mass loss by one of the stars (Blaauw 
1961). The latter is expected to occur in a supernova in 
which the exploding star loses a considerable fraction of its 
mass leaving behind a neutron star. The orbital periods of 
PSR B182Q-11 and PSR B2303-I-46 prohibit the observed 
eccentricity from being primordial; these binaries have ex- 
perienced a phase of mass transfer which circularized the 
orbit. The orbital eccentricities of these binaries are there- 
fore induced by a supernova which may have produced the 
currently observed pulsar. 

If the companions to PSR B1820-11 or PSR B2303-I-46 
are white dwarfs they must have been formed before the 
neutron stars. This possibility was never studied in detail 
since the progenitor of a white dwarf is expected to live 
longer than the progenitor of a neutron star. However, it 
was noticed by Tutukov & Yungelson (1993) and by Porte- 
gies Zwart & Verbunt (1996) that a reversal of evolution can 
be accomplished if both stars are similar in mass and if the 
masses of both stars are only slightly smaller than the mass 
limit for forming a neutron star. If this is the case, the sec- 
ondary star may still gain enough mass in a phase of mass 
transfer to pass the limit for forming a neutron star. For 
this to happen mass transfer should proceed rather conser- 
vatively, which will be the case since the initial mass ratio is 
close to unity. In the second, unstable, phase of mass trans- 
fer the white dwarf (original primary star) will spiral-in into 
the envelope of the secondary as this star ascends the gi- 
ant branch. This phase causes the orbital period to decrease 
dramatically by using orbital energy to carry the common 
envelope to infinity. If the white dwarf and the core of the 
giant stay detached after the common envelope is ejected, a 
close binary consisting of a white dwarf and a helium star 
remains. 

If the helium star is massive enough, the binary experi- 
ences a supernova and a young radio pulsar is formed. If the 
system survives the supernova the companion of the radio 
pulsar will be a white dwarf in an eccentric orbit. The mass 
of the white dwarf will be ^ LIMqIu such a case, because 
it originated from a star that was rather massive. 



2.3 Binary population synthesis 

We use the population synthesis program for binary stars 
SeBa^ (Portegies Zwart & Verbunt 1996) for evolving a 
million binaries with a primary mass between 8 Mq and 
100 Mq . We used model B from Portegies Zwart & Yun- 
gelson (1998; henceforth PZY98), which satisfactorily re- 
produces the properties of observed high-mass binary pul- 
sars (with neutron star companions). For a single star we 
adopt a minimum zero-age mass of 8 Mq for forming a neu- 
tron stars, which coincides with estimates based on observa- 
tions (Koester & Reimers 1996, see however Ritossa, Garcia- 
Berro, & Iben (1996) who show that a single star with a mass 



t The name SeBa is adopted from the Egyptian word for 'to 
teach', 'the door to knowledge' or '(multiple) star'. The exact 
meaning depends on the hieroglyphic spelling. 
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of 10 Mq may still evolve into a oxygen-neon-magnesium 
white dwarf). If a star is stripped of its envelope in an early 
phase this lower limit may become as large as 12 M0^. Co- 
incidence of the lower limits of the initial primary mass for 
our simulation with the minimum mass for the progenitors 
of neutron stars results in an underestimate in the forma- 
tion of binaries in which the secondary star experiences a 
supernova. By comparing our results with those of Portegies 
Zwart & Verbunt (1996, see their Tab. 4), who take lower 
mass primaries into account, we evaluate that this affects 
our results with 20%. 

While referring the reader to PZY98 for details, we 
list here the most crucial sissumptions. The masses of the 
primaries obey a power-law with exponent 2.5 (Salpeter 
= 2.35). The initial mass function is normalized to the 
current Galactic star- formation rate of 4MQyr~^ (van den 
Hoek & de Jong 1997). We assume that all stars are born in 
binaries with a semi-major axis up to a = lO'^ Rq. The dis- 
tribution was chosen to be flat in log a. Roche-lobe contact 
at zero age sets the minimum to the distribution of initial 
separations. The mass of the secondary is selected to be be- 
tween 0.1 Mq and the mass of the primary from a flat mass 
ratio distribution. The eccentricity for each binary is taken 
from the thermal distribution between zero and unity. 

Neutron stars receive a kick upon birth. The velocity 
of the kick is taken randomly from the distribution function 
proposed by Hartman (1997) and in a random direction. 
Black holes (from stars initially more massive than 40 M©) 
receive a kick velocity which is scaled by Mns/Mbh- 

From computations presented in PZY98 one may infer 
that our conclusions concerning birth rates and Porb — e dis- 
tributions for various combinations of neutron star and com- 
panion are robust with respect to reasonable variations in 
the most crucial model parameters and initial conditions, 
such as the kick velocity distribution and common envelope 
parameter, initial mass function and initial distribution in 
orbital periods. The criteria for the stability of mass trans- 
fer and the adopted rate of mass lost via a stellar wind and 
the amount of momentum lost per unit mass, however, may 
affect the details of our calculations rather signiflcantly. 

From the output of the computer simulations we select 
binaries with at least one neutron star. Within these con- 
straints three orbital-period ranges are considered: Porb > 
10, 100, and 1000 days. Table |l| gives the model birth rates 
for the binaries in various groups depending on the nature 
and the mass of the companion of the neutron star. Figures 
1 and 2 give the probability distributions for these binaries 
in the orbital period-eccentricity plane. 

The birth rates in Table |l| reflect the initial distribu- 
tions of the binary parameters and the complicated evolu- 
tionary history via several phases of mass transfer and at 
least one supernova. It is therefore not surprising that the 
most common companion to a pulsar is a rather massive 
( ^ 5M0) main-sequence star in a relatively wide binary 



i In SeBa the lower limit on the progenitor mass for forming a 
neutron star in a binary depends on the evolutionary stage at 
which the star loses its hydrogen envelope due to the interaction 
with its companion. This depends on the orbital separation and 
the initial mass ratio. The lower mass limit for the formation of a 
neutron star in a binary may thus be higher than for single stars. 



Table 1. Results of the model computations for binaries with 
at least one radio pulsar. The first column identifies the various 
binary components. Cols. (2) and (3) give the mass limits for the 
companion of the pulsar (black holes have a mass > 2M0). Col. 

(4) gives the total Galactic birth rate (BR) of such binaries. Cols. 

(5) , (6) and (7) give the birth rate of binaries with an orbital 
period larger than 10, 100, and 1000 days, respectively, in the 
units of 10~^ yr~^ 
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(10 ^ Porb ^ 1000 days. Fig. 1). As argued above, the pres- 
ence of such a massive companion in PSR B 1820-11 may 
be excluded by the absence of an optical counterpart and 
the unlikely inclination of the orbit; hiding a companion of 
M > 3Mq demands cosi > 0.95 (Lyne & McKenna 1989). 
The observed radio emission also makes it unlikely that the 
mass of the companion to PSR B1820-11 is ^ 5Mq. Such 
a companion star with a mass loss in the stellar wind of 
^ 10~® Mq yr"^ (which is chosen on the high side to ac- 
count for the enhanced wind mass loss for Be stars) would 
easily shield the radiation of the radio pulsar along its entire 
orbit (see Fig.|^). 

In the range of orbital periods Porb > 100, days the 
birth rates of neutron stars accompanied either by another 
neutron star (ns), a white dwarf (wd), a low mass (Af 
1.4 Mq) main-sequence star (ms) or a black hole (bh) are 
comparable (Table|^. The relatively large population of {bh, 
ns) systems and their distribution in orbital periods is a 
result of the evolutionary history of massive stars which is 
dominated by stellar wind mass loss rather than Roche lobe 
overflow and by the lower kick velocities imparted to black 
holes. However, the mass function of PSR B1820-11 makes 
a black hole with mass exceeding several Mq as a companion 
unlikely, and for PSR B2303+46 such a massive companion 
is excluded by the measured total mass. 

Figure |^ compares the probability distributions for bi- 
naries with a radio pulsar which is accompanied by a white 
dwarf (upper panel) and an older neutron star (lower panel) 
in the orbital period - eccentricity diagram. 

Our model calculations reveal that the binaries in which 
a young pulsar is accompanied by a massive white dwarf are 
generally rather eccentric (e ~ 0.6) and that their orbital 
periods are distributed with more or less equal probabil- 
ity in log Porb between a few hours and about one hundred 
days. In Fig. § PSR B1820-11 appears at the edge of the 
area with highest probability, while the orbital parameters 
of PSR B2303-I-46 fall nicely in the densely populated area 
of the Porb — e diagram. This suggests, together with the 
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Figure 1. Probability distribution in the orbital period — eccen- 
tricity plane for (ns, ms) binaries with a primary mass between 
O.7M0 and 1.4M0(upper panel), 1.4Mq and 5 (middle) and 
between 5 M-q and 10 Mq (lower panel, see Tab. n for the birth 
rates to which the panels are normalized). The solid line in the 
lower panel is plotted using Eq. (|^) for a 5 Mq star with a wind 
mass loss rate of lO~^M0yr~l. A binary is visible over its en- 
tire orbit to the right of both lines. To the left of both lines, the 
binary stays invisible throughout its entire orbit. In between the 
left and the right part of the solid lines the binary is visible for 
part of its orbit (near apocenter). The dashed line is plotted using 
a 10 Mq primary with a mass loss rate of lO^'^MQyr"^. The * 
symbol indicates the position of PSR B1820-11, and two bullets 
(lower panel) indicate the positions of PSR B1259-63 (right) and 
J0045-7319 (left). 



rather high birth rate of such binaries, that PSR B2303+46 
may weU contain a white dwarf as a companion star. Re- 
cently van Kerkwijk & Kulkarni (1999) found evidence that 
PSR B2303-I-46 is accompanied by a white dwarf with a 
mass of about 1.2 Mq. 

PSR B 1820-11 fits rather iU in either of the Porb - e di- 
agrams for (ns, ns), {wd, ns), and low mass or moderate 
mass (ns, ms) binaries. The relative birthrates for a popu- 
lation of binaries with long orbital periods increases by in- 
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Figure 2. Probability distributions for binaries which contain a 
young neutron star and a white dwarf (top panel) or a neutron 
star (lower panel) as a companion star. A figure for {bh, ns) is 
not provided because of the small number of data points. The 
* symbols identify the positions of PSR B1820-11 (right) and 
PSR B2303-I-46 (left). Dots in the lower panel give the positions 
of the known recycled pulsars with a neutron star companion. 



creasing the efficiency of the deposition of orbital energy into 
common envelopes during the spiral-in phase (see PZY98). 

However, this would require us to relax the decisive 
criterion why we chose this specific set of model parame- 
ters; the coincidence of having the other high mass binary 
pulsars (PSR J1518-H49, B1534-hl2, and B1913-hl6) in the 
most densely populated areas in the Porb — e diagram for 
(ns, ns) binaries. There is no evidence for a strong selection 
against the discovery of a young pulsar which is accompa- 
nied by a dead pulsar or by a massive white dwarf with an 
orbital period ^ 10 days (Johnston & Kulkarni 1991). The 
higher birth rate of binaries in which the radio pulsar is ac- 
companied by a low mass main-sequence star and the better 
coincidence in the Porb ~ e diagram makes this possibility 
very attractive. 

If the companion of PSR B1820-11 is indeed a low- 
or moderate mass main-sequence star it may be a precur- 
sor of a system similar to GX l-|-4 = V2116 Oph, which 
shows the X-ray features of an accreting neutron star and 
which is identified with a symbiotic star (Davidsen, Malina, 
& Bowyer 1977). The orbital period of GX 1-1-4 is > 100 
days (Chakrabarty & Roche 1997). The low birth rate of 
such binaries (Table |l|) together with the short life time of 
a red giant explains the paucity of such X-ray binaries. 
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3 CONCLUSION 

We have studied the nature of the companions of young radio 
pulsars in eccentric orbits. Our calculations reveal that the 
possibility that a young pulsar is accompanied by a white 
dwarf of mass ^ I.IM0 cannot be neglected. The birth rate 
of such binaries is higher than the birth rate of of binaries in 
which a young pulsar is accompanied by a dead pulsar (i.e.: 
an old neutron star). 

We argue therefore that PSR B 1820-11 may be accom- 
panied by a massive white dwarf, however the probability 
that its companion is a main-sequence star with a mass 
^ 5Mq (or even another neutron star) is similar to the prob- 
ability of being a white dwarf. 

Our calculations confirm the observation of van Kerk- 
wijk & Kulkarni (1999) in the sense that it is likely that 
PSR B2303-I-46 is accompanied by a white dwarf with a 
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